Localized many particle Majorana fermions with vanishing time reversal symmetry breaking in 

double quantum dots 
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We introduce the concept of spinful many-particle Majorana fermions. These emergent particles are nonlocal 
in even operator products, but local in odd operator products, thereby preserving their local statistics. We 
consider a superconductor - double quantum dot system where these modes can arise in the presence of vanishing 
Zeeman splitting. We find a reverse Mott-insulator transition, where the even parity band crosses the odd parity 
band. Above this transition, Majorana operators move the system between the odd parity ground state, associated 
with elastic co-tunneling, and the even parity ground state, associated with crossed Andreev reflection. These 
Majorana modes are described in terms of one, three and five operator products. Parity conservation results in 
a An periodic supercurrent in the even state and no supercurrent in the odd state. The necessity of only small 
magnetic fields and the broad parameter space where Majorana fermions appear simplifies the experimental 
realization greatly, paving the way towards the detection of the Majorana fermion. 



The prediction for the existence of Majorana fermions [|lj 
has attracted enormous attention in condensed matter physics 
in recent years Allured by the possibility of con- 

structing topological qubits for quantum computation 13,0], 
a plethora of schemes promising the positive identification 
of Majorana fermions has emerged fi- 12|. Superconductors 
provide intrinsic electron-hole coupling resulting in charge- 
less quasiparticles, and most proposals to realize Majorana 
fermions are based on this principle. Standard s-wave super- 
conductors do not support a zero energy Majorana mode. In 
the presence of strong Rashba spin-orbit coupling, however, 
an effective spinless p-wave superconductor can arise. This 
idea was first recognized by Fu and Kane, who considered 
the interface of a superconductor and a topological insulator 
Istl- A topological state harnessing Majorana modes can also 
arise using semiconducting nanowires in the presence of Zee- 
man and Rashba fields laS]- Experimentally, sup ercurrents 
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1 13], Fraunhofer patterns and Shapiro steps 111411 . SQUIDs 
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and zero bias conductance peaks Ill6[l have been ob- 
served in topological insulator systems, which together with 
the zero bias conductance peaks in nanowire systems lIlTIl pro- 
vide prospects for the observation of the Majorana fermion, 
but to date, no conclusive evidence has been observed. 



Quantum dots can also be used to realize localized Ma- 
jorana modes 111 ll llSl 11911 . Quantum dot proposals include 



the use of Rashba coupling 111 111 , but the presence of an 
anisotropic magnetic field can also mimic this effect lfl9ll . 
Although these proposals are already experimentally quite 
promising, the strong Rashba or anisotropic magnetic fields 
required makes them very challenging and limits material 
flexibility. For example, anisotropic magnetic fields only re- 
sult in spinless localized Majorana modes in fields Ez ^ A, 
with A the induced superconducting gap. Nonetheless, here 
we will show that in the presence of small anisotropic mag- 



FIG. 1 : (a) Schematic representation of the device considered. Two 
superconductors (with phase difference (p-) are connected via a dou- 
ble quantum dot with on-site Coulomb repulsion U and in the pres- 
ence of a magnetic field Ez.± ~ fcsT. A nanomagnet introduces a 
localized Ez^\\, which rotates the field near one of the two dots, and 
we define 6 to be the relative angle between the local fields, (b) Tun- 
ing to the degeneracy of the ground state of the double quantum dot 
system. The color scale represents the energy gap between the lowest 
even and odd energy eigenstates. Eg = (ei— e2)/A. The parameters 
are U = 3|A| — —3t. The degenerate ground states support many- 
particle Majorana fermions with local statistics, for all 9 0. The 
system is symmetric around = tt, but not around 6 = n/2, since 
the presence of regular Andreev reflection distinquishes the quantum 
dots having parallel and antiparallel spin axes. 



netic fields, \B\ ^ ksT, a new arena emerges: the concept 
of spinful many-particle Majorana fermions. These Majorana 
fermions are still localized in their odd operator products, and 
thereby preserve the local statistics of spinless proposals 12111 . 

We consider an s-wave superconductor - double quantum 
dot system as depicted in Fig. IT], where crossed Andreev re- 
flection (CAR) dominates elastic co-tunneling (EC), a regime 
that is readily achieved 122-261. Double electron occupancy 
in a quantum dot, as we will show, inhibits, but does not pre- 
vent the appearance of Majorana fermions. On-site Coulomb 
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repulsion U can be used to tune the system from the dou- 
ble occupancy regime, which we define as [/ = 0, towards 
the single electron regime, U ^ oo. Upon increasing the 
Coulomb repulsion, we find a clear phase transition where it 
becomes possible to tune the even and odd parity states to be- 
come degenerate, see Fig.||2l. This is reminiscent of a Mott 
metal-insulator transition, however in reverse, where the on- 
site Coulomb repulsion drives the system toward having a de- 
generate ground state. Above this transition, we can always 
tune towards the degeneracy point using a second supercon- 
ductor with a phase difference Many-particle Majorana 
fermions appear right at the crossing and are described by one, 
three and five operator products. The even operator products 
are nonlocal, while the odd operator products are local for all 
angles 9, the angle between the local fields on the dots. 

The requirements of only small fields and the use of quan- 
tum dots opens large material flexibillity. For example, the 
system could be based on III-V heterostructures 1271], but also 
on materials such as silicon and diamond with large fractions 
of non-magnetic nuclei which increase spin coherence times 
lEill . In order to elucidate the novel many -particle physics that 
emerges, we will take all results in the limit Ez 0, valid 
as T — > 0, meaning we have spinful, spin rotation symmetric, 
ground states. 

The superconductor - double quantum dot system with on- 
site Coulomb repulsion in the presence of an anisotropic mag- 
netic field, Fig.ITl, is described by the Hamiltionian 



H — Hs + Hd + H' 
The superconducting part, Hs, is 
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(1) 



Here, a denotes the spin that is not a. The double quantum 
dot is described by, Hd, 

2 

with U the onsite Coulomb repulsion, and tj the onsite energy 
of the j*^ dot. The third term introduces the tunneling and is 
described by, Ht, 



(4) 



The overlap integral Fy is between the i quasiparticle in the 
superconducting lead, and the j"^ dot. Finally, 



of Ref.|29fl, to obtain the effective Hamiltonian of the double 
dot-superconductor system. 

Before presenting the final effective theory, we first con- 
sider the anisotropic magnetic field, which defines the spin 
axes. The angle 9, the angle between the two local magnetic 
fields, modify EC and CAR as follows: 

tc\ „C2,cr tcos{9/2)c{ ,^C2,a + s'm{9 /2)c{^^C2,g- 
Aci,crC2,ff crAsin(6'/2)ci,CTC2,CT + A cos(6'/2)ci,ctC2,s-, 

(6) 

where cr, a = ±, and t, A are the effective EC hopping and 
CAR Cooper pairing amplitudes 1291]. The final effective 
Hamiltonian of the dot-superconductor system is 



^(i cos(0/2)cl^,C2,a + atsin(0/2)cl_^C2,a 



(^) {a sin(0/2)ct4,, + cos{9 l2)clA,s) 



~Ae* 2 cos 



+Ae'^ cos(^) Y AA,s + 



(7) 



where <\)± = 0i ± 02 is the sum (difference) between the 
phases of the two superconductors. From hereon we will as- 
sume that the onsite energy of the two dots has been tuned to 
the chemical potential of the superconductors, which we de- 
fine as our zero of energy (ei = £2 = 0). A discussion of the 
effects of the onsite energies deviating from this 'sweet spot' 
has been presented elsewhere [fl9ll . 

The Hamiltonian Eq. ]|7] cannot be decomposed into to- 
tal spin sectors, as the anisotropic magnetic field mixes these. 
However, the fermion parity is a conserved quantity. The to- 
tal Hilbert space has dimension (2"'')"j /2 = 8 for each par- 
ity sector (i.e. even and odd), which makes exact diagonal- 
ization particularly straight-forward. In the occupation rep- 
resentation, we can define an occupation basis |1,2), where 
the numbers correspond to the two dots, and we will use ar- 
rows to denote the spin, and then we can construct a basis 
of the sixteen possible configurations. In the [/ — > 00 limit, 
the total Fock space is restricted to nine possible states, 5 for 
even parity, and 4 for the odd, which can be easily solved 
to give the eigenstates that determine the ground state, in the 
even and odd sectors. We solve for finite Zeeman splitting, 
and obtain in the odd parity sectors that the eigenstates are 
todd = ±i ± Ez- We take the limit Ez — > in the low- 
est energy odd parity state to obtain the limit of vanishingly 
small magnetic fields of the ground state even and odd parity 
wavefunctions 



(5) 



The Hamiltonian is quadratic in the leads, and so we can inte- 
grate out the superconducting leads, following the procedure 



±y2e^|0,0) +cos(0/2)(| t,!) + I ;,t» 
+ sin((?/2)(|t,t)-U,;)) 



(8) 
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FIG. 2: Phase diagram showing the onset of a degenerate ground 
state energy gap as a function of on-site Coulomb repulsion, for 
|A| — \t\. Upon increasing the on-site Coulomb repulsion U, a re- 
verse Mott-insulator transition occurs where the lowest energy state 
changes from even to odd parity. At the crossing point, Majorana 
operators move the system between the degenerate ground states. 



with corresponding eigenvalue 
and 



TV2Acos{(f)^/2), 



=in(0/4)(|t,O) + |O,t)) 

=os(0/4)(|;,o)-|o,;))" 



(9) 



with eigenvalue eodd = ~t (taking t > 0). A degenerate 
ground state is obtained when feuen — ^odd- Crucially, for 
any t < \/2\A\, there is always a (/)_ which can be chosen, 
depending on the value of 9, such that a degenerate ground 
state can be obtained. 

We emphasize that as we are using the occupation number 
basis, the ground state is given by the lowest energy eigen- 
state, and not the zero energy eigenstate as in the more famil- 
iar Bogolubov de-Gennes theory. 

The two degenerate ground states have even and odd parity, 
respectively, and so are protected from hybridizing when the 
total system conserves particle number parity. When there are 
two degenerate ground states, a pair of Majorana operators 
(71 , 72) can be constructed which transform the two ground 
states into each other, such that 71^'odd = '^even, for exam- 
ple. These Majorana operators are given by 
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where x = (sin(6'/4) — cos(0/4)) ^, and 
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FIG. 3: Josephson supercurrent of the even parity sector at T = 
0.01 1 A| for the limiting cases U = Q and cxd and intermediate 
U = 20|A|, with |A| = \t\ and 6 = 7r/2. When C/ = 0, reg- 
ular Andreev reflection dominates and the Josephson current is 27r 
periodic and no Majorana fermions exist at all. When U — > 00, reg- 
ular Andreev reflection is completely forbidden, there zero energy 
states are local odd operator Majorana fermions, and crossed An- 
dreev reflection results in a 47r periodic Josephson effect. For finite 
U , hybridization between the even parity eigenstates leads to a sharp 
transition at <^ = tt. For large U, Zener tunneling F can restore the 
A-K periodicity, but Majorana modes are always present. Jo = eA//i 
is the maximum supercurrent at [/ = 0. The odd parity supercurrent 
is always 27r periodic, and small, vanishing completely in the limit 
U ^ 00. 



and n, 



(11) 



' (6*) has a similar form with site 
indices 1 , 2 interchanged. Eq. ifTTl has the form of a usual 
Majorana operator 112 ill , except that the phase dependence is 
the total phase of the two superconductors divided by four, or 
equivalently half the average phase of the two superconduc- 
tors. 

The form of three-operator Majorana fermion in Eq. [ [TOl 
is nonlocal in general. However, it is interesting to note that 
the relative phase of the Majorana components, which is re- 
sponsible in general for their non-Abelian braiding statistics, 
is still localized, as the nonlocality of the Majoranas as defined 
above is restricted to number operators only. It is assumed, in 
the above Majorana expressions, that the system is projected 
into the singly occupied Pock space, since J7 — )- 00. 

When U is finite, a phase transition occurs at a critical 
value of the on-site repulsion. For Coulomb repulsions be- 
low this critical point, the ground state is non-degenerate, with 
the lowest energy even parity state always having lower en- 
ergy than the lowest energy odd parity state. In Fig.H] we 
have plotted the excitation energy of the first excited state 
Eg = (£2— ei)/|A| as a function of onsite Coulomb repulsion, 
and relative spin angle 6. The phase transition corresponds to 
the vanishing of the excitation gap at a critical value of on-site 
Coulomb repulsion. This critical value varies as a function of 
angle, but for aU angles it is quite small, being of the same 
order as the effective pairing | A | . 

For finite on-site Coulomb repulsion, the structure of the 
ground state, and Majorana states are much more complicated 
than in the infinite U case, though qualitatively do not differ 
The odd parity ground state develops a finite weighting on the 
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terms of the form | tii t)' together with the three equivalent 
combinations of this, whilst the even parity ground state de- 
velops a finite weighting on doubly occupied dots (| '11,0) 
and |0, ti))' together with the four fermion, double occupied 
quantum dot pair | The Majorana fermions, how- 

ever, retain the same phase locality condition. The Majo- 
rana fermions in this case acquire five fermion operator prod- 
ucts, together with single and three fermion operator products. 
Generically, the Majorana has the form 




E 

aa'jj' 



(12) 

where 6^ ^ = <,„m = = 0, and when (/> = tt/2, 

all operators where j,j' = 2 have coefficient zero. The ap- 
pearance of three and five operator products are a clear gen- 
eralization of the Majorana fermion concept, which is usually 
based on single operator products |23 21 1. 

The regime where Majorana fermions appear can be found 
by analyzing the Josepshon supercurrent through the double 
dot system. The Josephson current is calculated from the 
derivative of the free ener gy w ith respect to the superconduct- 
ing phase difference 13011 . In contrast to topological su- 
perconductor systems, the odd parity ground states disperse 
only very weakly with superconducting phase difference, and 
are strictly 27r periodic. At J7 — > oo, they are completely flat, 
as CAR cannot possibly excite the odd parity ground state. 
The absence of supercurrent in the odd parity state will be a 
strong signature of parity consen'ation. In Fig.fSl we have 
plotted the Josephson current for the even parity system as a 
function of superconducting phase difference for the U oo 
case, where Majorana fermions are present (i.e. CAR splits 
the Cooper pairs and a gapless super-current flows between 
the superconductors) versus the case of regular Andreev re- 
flection where the final term in Eq. ||7] dominates, as in the 
U = case. The periodicity of the two cases is clearly visi- 
ble. We have also plotted an intermediate case, U ~ 20|A|, 
which is 27r periodic in the adiabatic limit. We note that 
while in a p-wave superconductor the boundstates are pro- 
tected by fermion parity, the dispersive bound states in the 
double quantum dot system have both even parity, and scat- 
tering between the branches is possible even when the total 
parity is conserved. Interestingly, the matrix element coupling 
the branches is a two Majorana operator product. The anoma- 
lous current phase relationship, however, can still be observed 
in non-equilibrium measurements [31] or using dc SQUIDs 
llllll. 

At finite U, the bonding and anti-bonding even parity states 
hybridize as (/>_ drives them toward a degeneracy point, nat- 
urally resulting in a signal where the bonding state tunnels 
into the anti-bonding, and one where it does not, as shown in 
Fig- El- A superposition of a 27r periodic signal and a Air pe- 
riodic signal is expected in this case, since by means of Zener 
tunnelingT the quasiparticle can overcome the hybridization 
gap I34II35II . This transition from 4tt to In periodicity is not 
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FIG. 4: Energy eigenstates of the double dot system for [/=20|A|, 
|A| = |t| and 9 = n/2. The +, — signs indicate the even par- 
ity ground state bonding and antibonding state, determined by the 
relative sign between the empty and two-fermion occupations, for 
example *± = a|0000) ± 6|1010) + .... At infinite U, the two 
states are orthogonal, but at finite U the two develop an anticross- 
ing. The arrow indicates the evolution of the ground state that gives 
a 4-K periodic Josephson effect, which is still possible for small hy- 
bridization gaps by means of Zener tunneling F. As the bonding and 
anti-bonding states become increasingly hybridized, the 2n periodic 
Josephson effect becomes dominant. We emphasize that as we are 
using the occupation number basis, the ground state is given by the 
lowest energy eigenstate, and not the zero energy eigenstate as in the 
more familiar Bogolubov de-Gennes theory. 



a transition where Majorana modes disappear, but a transition 
where the five operator products vanish. Also, finite U results 
in a weakly dispersing odd parity ground state, clearly visi- 
ble in Fig.lH, and introduces a supercurrent in the odd parity 
state. 

In conclusion, we have presented an experimentally fea- 
sible scenario whereby many particle Majorana fermions in 
double quantum dots may be realized. We have shown that 
a large Zeeman splitting is not necessary to realize nonlocal 
Majorana fermions in double quantum dots, due to the onsite 
Coulomb repulsion, and in fact a small Zeeman field is suffi- 
cient. Furthermore, we have shown that by connecting two out 
of phase superconductors to both quantum dots, the parame- 
ter space can be tuned so as to obtain a degenerate ground 
state protected by parity, with Majorana operators connect- 
ing the ground state. The Majorana operators are local in the 
odd operators, which determine their statistics, for any angle 
6. The effect of a finite, rather than infinite, on-site interac- 
tion does not affect the locality conditions of the Majoranas, 
nor their parity. We expect that the proposed device will ex- 
pand Majorana physics even further, not in the least due to 
its spinful ground state, and the fact that the corresponding 
Majorana fermions are constructed from one, three, and five 
fermion creation/annihilation operator products, together with 
their Hermitian conjugates. 
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